Possible evidence for the mass shift of 77' meson 
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Topological charge distributions of SU(3) gauge theory at finite 

Q_i' temperature are calculated on the lattice with high statistics in a 

^SZ3 ■ manner of free from the uncertainty of /3-f unction. Clear tempera- 



ture dependence of the topological charge distribution is obtained. 
Combining the present result with the recent study about the re- 
lation between the masses of r\ meson system and the gauge field 
topology, we obtain m v ,{T)/m v ,{0) = 0.86 + 0.02 and 0.69 + 0.01 for 
T = 0.75T C and T = 0.93T C , respectively. The result shows the first 
clear evidence of the mass shift of rf meson at finite temperature. 



1 Introduction 

It is a common knowledge that at sufficiently high temperature, 
the hadronic system may undergo a phase transition to the QGP 
(quark-gluon-plasma) phase. What is a signature or a precursor of 
the QGP is a long standing question. Although a lot of possibilities 
have been posted so far, there is no definite answer yet Jl|. Among 
them the mass shift of hadrons might be an promising candidate for 
the precursor of QGP if this phenomena is in the case. 0, |3|, [|, |5|] 

This is because the mass shift, itself, can indicate that the system 
is very close to the critical region without a precise quantitative 
argument. 

In this paper, we focus on the neutral pseudo-scalar meson system 
whose mass spectrum at zero temperature can be well explained by 
the linear sigma model incorporated with the axial U(l) anomaly 
||. It is also well known that the axial U(l) anomaly is caused by 
the topological fluctuations of the system, 

df+j^x) oc Q(x), (1) 

where j^ is the flavor singlet axial vector current and Q(x) is the 
topological charge density. 

Recently the relation between the singlet-octet mass difference of 
the neutral pseudo-scalar meson system and the topology has been 
analyzed in detail and the strong correlations of these two has been 
suggested 0. 

From these results mentioned above, if there exist the tempera- 
ture dependence in the topological sector of the theory, it might be 



reflected by the possible change of the mass spectrum. 

The organization of this paper is as follows. In the next sec- 
tion 2, we briefly mention our model of the neutral pseudo-scalar 
meson system as well as it's relation with the gauge field topology. 
The lattice results of the topological charge distributions at finite 
temperature is given in section 3. In the section 4, we will discuss 
the possible shift of the masses. The section 5 is devoted to the 
conclusion. 

2 Model of the neutral pseudo-scalar meson sys- 
tem 

For the model of the neutral pseudo-scalar meson system, we adopt 
the following mass matrix incorporated with the axial U(l) anomaly 
which was introduced by Pisarski and Wilczek || and was able to 
reproduce the real mass spectrum at zero temperature. 
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where M 2 is the (mass) 2 matrix for the ttq — tj — t]' system, rrii is the 
quark mass for i-quark (i=up, down and strange), fj is the decay 
constant for j-meson (j=7i"o, r\ and rj'), v is the strength of the chiral 
condensate (— v =< uu >=< dd >=< ss >) and K is the free 
parameter representing the effects of the axial U(l) anomaly. 

In order to solve this equation at finite temperature, the temper- 
ature dependence of all parameters in the above equation should be 
given. Pisarski and Wilczek have considered only the K effect of 
the system. [|| Although there have been several model works which 
suggest the possible change of other parameters like /„. and v at 
finite temperature, the recent lattice calculation indicates that the 
temperature dependence of these quantities seems very small up to 
very close to the critical temperature 0. 

So in this paper we also restrict ourselves on the K effect of the 
system. 

Let us first consider the case that the flavor SU(3) symmetry is 
unbroken, i.e., m u = m^ = m s . In this case, K is nothing but the 
mass difference mo between singlet and octet meson which is directly 
related the topological charge by the result of ref . . Then once we 
know the topological charge distribution p(Q), we can estimate K 
as. 

K = m = JdQp(Q)m (Q), (3) 

where 

p(Q)=p(Q)/JdQp(Q) (4) 

is the normalized topological charge distribution. 

In eq.2, it is also clear that K represents the strength of the axial 



U(l) anomaly which is connected with the gauge field topology via 
eq.l. So it is rather plausible to assume that the K is independent of 
the flavor SU(3) symmetry. Under this assumption, above equation 
becomes valid even in the flavor symmetry broken phase. 

Next we assume the singlet-octet mass difference m is purely 
due to the gauge field topology, i.e., all temperature effect to m 
should be included through the change of the topological charge 
distribution p(Q). Then we can use eq.3 even at finite temperature. 

3 Topological charge distribution at T^O 

When we evaluate eq.3 at finite temperature on the lattice, there 
are two restrictions at this stage; 

• mo(Q) was obtained at only one (3{— 6/g 2 ) 0, 

• unknown non-perturbative behavior of /3-function. 

In order to get rid off these restrictions, we have decided to make 
simulation at fixed f3{= 6/g 2 ). In this case, T^O can be realized by 
changing N t (= temporal extent of the lattice) through T = 1/N t a. 
( a is the lattice spacing for all directions. ) 

The simulation is done on 16 3 x 16, 8 and 6 lattice in SU(3) 
gauge at (3 = 5.89. These lattices correspond to the temperature 
T ~ 0, 0.75T C and 0.93T C , respectively, by the fact that the critical 
(3 on 16 3 x 6 lattice is (3 C ~ 5.894 || and by using a -1 vs. (3 data 
by QCDTARO Collaboration |II|. For the SU(3) simulation code 
in this analysis, we use QCDMPI M. 



After the thermalization of 1000 pseudo-heat-bath sweeps, the 
measurement is done every 100 sweeps. Number of configurations 
used are 3500, 1500 and 3200 for T ~ 0, T = 0.75T C and T = 0.93T C , 
respectively. To extract the topological charge Q, we adopt the 
cooling method [ |TT|j . 

The result of p(Q) is shown in Fig.l. In this figure, clear tempera- 
ture dependence of p(Q) is seen. The distribution becomes narrower 
as the temperature gets close to the critical temperature. This fea- 
ture can be understood as the partial restoration of the axial U(l) 
symmetry at finite temperature. 

The change of the topological charge distribution p(Q) at finite 
temperature should be treated carefully. If we assume the distribu- 
tion is Gaussian, p(Q) oc exp(— aQ 2 ), the constancy of the topolog- 
ical susceptibility in the confinement phase |T2j requires, 



a oc N- 1 (5) 

at fixed p. 

In this sense, the distribution becomes narrower as the tempera- 
ture increases (N t decreases). So this factor (eq.|j) should be taken 
into account for the quantitative study in the next section. 

4 Temperature dependence of K and m v 

For the parameterization of m (Q), we adopt two different fits of 
the result in ref.@. One is the least square linear fit (fit A); 

mo(|Q|)=Ci|Q| + C2, (6) 



with c\ = 0.08 and C2 = 0.20 in MeV units 0. It is surprising that 
this linear fit seems very well as you can see in ref . . In this fit C2 
is finite. However it is natural for C2 to be zero at |Q|=0 as far as 
the flavor SU(3) symmetry is unbroken. So we make another fit (fit 

B); 

m (\Q\) = c[\Q\ (7) 

with c[ = 0.16, such that the average mass splitting is almost equal 
for these two fits. 

Since f3 in this simulation is different from (3 in ||, we can not fix 
the absolute scale without the knowledge of the /3-function of SU(3) 
gauge theory. So we concentrate on the mass ratio normalized to the 
zero temperature value which is independent of this (3 discrepancy 
as far as the fits A and B are concerned. 

When we calculate K(T) in eq.|3|, N t effect in eq.^| should be 
eliminated. To do this, in this paper, we normalize K(T) by its 
Gaussian value, 

js(rp\ -^ {■*■ J lattice / \ 

K K T ) = jTrn • ( 8 ) 

-**■ V Jgaussian 

Combining mo(Q) and eq.3, we then obtain the temperature de- 
pendence of K parameter. The result is shown in Fig.2. We can see 
the dramatic decrease of this quantity as the temperature gets close 
to T c . 

Finally we solve the eigenvalue problem of eq.2 using this K(T). 
Other parameters are taken from ref. 0. The result for m v (T) and 
m^iT) are shown in Fig. 3. In these figures, only statistical error 
coming from the topological charge measurement has been taken 
into account. Although systematic error coming from above 2 fits 



(fit A and fit B) is not so small, mass decrease near T c is seen for 
both r) and 77'. (m n does not show any noticeable change in this 
calculation as in ref. |§.) Especially for rj', clear mass shift can be 
seen (in the case of fit B), 

m ,l . = 0.86 ± 0.02(T = 0.75T C ) (9) 

77ty(T=0) V ' V ^ 

= 0.69 ± 0.01 (T = 0.93T C ) (10) 

This feature can also be understood as the direct reflection of the 
weakenness of the axial U(l) anomaly. From the present result, 
if we measure this temperature dependence in the experiment, it 
could be a possible precursor and/or thermometer of the QGP phase 
transition. 

5 Conclusion 

We have analyzed the temperature dependence of the topological 
charge distributions of the SU(3) gauge theory. Measurement has 
been performed on the lattice at fixed /3(= Q/g 2 ) which enable us 
to evaluate the temperature dependence of the physical quantity 
free from the uncertainty of /^-function. Although the systematic 
error coming from the uncertainty of mo(Q) is not so small, we 
find the clear mass shift of 7/ meson near the critical temperature. 
This is a first clear evidence of the mass shift of 7/ meson at finite 
temperature. 

There are several assumptions in this paper. Although these 
assumptions seem plausible at this stage, the detail investigation are 
highly expected in order to confirm the present result. Calculations 



with dynamical quark effect should be performed in future because 
the dynamical quarks may play a crucial role on this problem. 
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figure caption 

• Figure 1 

Topological charge distributions p(Q) at (a) T ~ 0, (b) T = 
0.75T C and (c) T = 0.93T C . 

• Figure 2 

K(T) normalized to the zero temperature value. • is for m (|Q|) 
c i\Q\ + C2 fitting and o is for m (|Q|) = c[\Q\ fitting. Errorbars 
represent only statistical errors. 

• Figure 3 

(a) 77' and (b) rj mass vs. temperature normalized to the 
zero temperature value using the mo(\Q\) result by Fukugita 
et al.0, where m is the r/ '-octet mass splitting. • is for 
m o(|Q|) = c i\Q\ + c 2 fitting and o is for mo(|Q|) = c[\Q\ fitting. 
Errorbars represent only statistical errors. 
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Figure 2 
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Figure 3(a 
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Figure 3(b 
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